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The effects of random magnetic fields are considered in an Ising spin-glass model 
defined in the limit of infinite-range interactions. The probability distribution for 
the random magnetic fields is a double Gaussian, which consists of two Gaussian 
distributions centered respectively, at +Hq and —Hq, presenting the same width a. 
It is argued that such a distribution is more appropriate for a theoretical description 
of real systems than its simpler particular two well-known limits, namely the single 
Gaussian distribution {a ^ ^o); and the bimodal one (cr = 0). The model is 
investigated by means of the replica method, and phase diagrams are obtained within 
the replica-symmetric solution. Critical frontiers exhibiting tricritical points occur 
for different values of o", with the possibility of two tricritical points along the same 
critical frontier. To our knowledge, it is the first time that such a behavior is verified 
for a spin-glass model in the presence of a continuous-distribution random field, 
which represents a typical situation of a real system. The stability of the replica- 
symmetric solution is analyzed, and the usual Almeida-Thouless instability is verified 
for low temperatures. It is verified that, the higher-temperature tricritical point 
always appears in the region of stability of the replica-symmetric solution; a condition 
involving the parameters Hq and o", for the occurrence of this tricritical point only, 
is obtained analytically. Some of our results are discussed in view of experimental 
measurements available in the literature. 
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I. INTRODUCTION 



Spin-glass systems y, 0, [s] continue to challenge many researchers in the area of 
magnetism. /^From the theoretical point of view, its simplest version defined in terms of 
Ising spin variables, the so-called Ising spin glass (ISG), represents one of the most fasci- 
nating problems in the physics of disordered magnets. The ISG mean-field solution, based 
on the infinite-range-interaction model, as proposed by Sherrington-Kirkpatrick (SK) j^, 
presents a quite nontrivial behavior. The correct low-temperature solution of the SK 
model is defined in terms of a continuous order-parameter function [3| (i.e., an infinite 
number of order parameters) associated with many low-energy states, a procedure which 
is usually denominated as replica-symmetry breaking (RSB). Furthermore, a transition in 
the presence of an external magnetic field, known as the Almeida-Thouless (AT) line js], 
is found in the solution of the SK model: such a line separates a low-temperature region, 
characterized by RSB, from a high-temperature one, where a simple one-parameter solu- 
tion, denominated as replica-symmetric (RS) solution, is stable. The validity of the results 
of the SK model for the description of more realistic systems, characterized by short-range- 
interactions, represents a very polemic question j^. Recent numerical simulations claim 
the absence of an AT line in the three-dimensional short-range ISG js], as well as along the 
non-mean-field region of a one-dimensional ISG characterized by long-range interactions 
[lo| . However, these results, obtained with rather small lattice-size simulations, do not 
rule out the possibility of a crossover to a different scenario at much larger lattice sizes, 
or also for smaller fields (and/or temperatures). One candidate for alternative theory 
to the SK model is the droplet model till, based on domain- wall renormalization-group 

nn 

arguments for spin glasses [12, Il3| . According to the droplet model, the low-temperature 
phase of any finite-dimensional short-range spin glass should be described in terms of 
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a single thermodynamic state (together, of course, with its corresponding time-reversed 
counterpart), i.e., essentially a RS-type of solution. Many important features of the ISG 
still deserve an appropriate understanding within the droplet-model scenario, and in par- 
ticular, the validity of this model becomes questionable for increasing dimensionalities, 
where one expects the existence of a finite upper critical dimension, above which the 
mean-field picture should prevail. 

Some diluted antiferromagnets, like FexZni_xF2, FexMgi_xCl2 and MnxZni_xF2, have 
jeen the object of extensive experimental research, due to their intriguing properties 
1^ . These systems are able to exhibit, within certain concentration ranges, random- 
field, spin-glass or both behaviors, and in particular, the compounds FexZni_xF2 and 
FexMgi_xCl2 are characterized by large crystal-field anisotropics, in such a way that they 
may be reasonably well-described in terms of Ising variables. Therefore, they are usually 
considered as good physical realizations of the random- field Ising model (RFIM), or also 
of an ISG. For the FexZni_xF2, one gets a RFIM-like behavior for x > 0.42, an ISG for 
X ~ 0.25, whereas for intermediate concentrations (0.25 < x < 0.42) one may observe 
both behaviors depending on the magnitude of the applied external magnetic field [R FIM 



(ISG) for small (large) magnetic fields], with a crossover between them [15|, ll6|, ll7|. In 
what concerns FexMgi_xCl2, one gets an ISG-like behavior for x < 0.55, whereas for 0.7 < 
X < 1.0 one has a typical RFIM with a first-order transition turning into a continuous 
one due to a change in the random fields Q, Even though a lot of experimental 
data is available for these systems, they still deserve an appropriate understanding, with 
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only a few theoretical models proposed for that purpose 
Within the numerical-simulation technique, one has tried to take into account the basic 
microscopic ingredients of such systems [20|, 
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231], whereas at the mean-field level. 



a joint study of both ISG and RFIM models has been shown to be a very promising 



approach [24 



25 



27|. 



In the present work we investigate the effects of random magnetic fields, following a 
continuous probability distribution, in an ISG model. The model is considered in the 
limit of infinite-range interactions, and the probability distribution for the random mag- 
netic fields is a double Gaussian, which consists of a sum of two independent Gaussian 
distributions. Such a distribution interpolates between the bimodal and the simple Gaus- 



4 



sian distributions, which are known to present distinct low-temperature critical behavior, 



within the mean-field limit 
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271]. It is argued that this distribution is more 



appropriate for a theoretical description of diluted antiferromagnets than the bimodal 
and Gaussian distributions. In particular, for given ranges of parameters, this distribu- 
tion presents two peaks, and satisfies the requirement of effective random fields varying 
in both sign and magnitude, which comes out naturally in the identification of the RFIM 
with diluted antiferromagnets in the presence of a uniform field j^, 29 1; this condition is 
not fulfilled by simple discrete probability distributions, e.g., the bimodal one, which is 
certainly very convenient from the theoretical point of view. Recently, the use a double- 



Gaussian distribution in the RFIM 



301] yielded interesting results, leading to a candidate 



model to describe the change of a first-order transition into a continuous one that occurs 
in FexMgi_xCl2 [l4, Q, 19 1. The use of this distribution in the study of the present 
model should be relevant for FexMgi_xCl2 with concentrations x < 0.55, where the ISG 
behavior shows up. In the next section we study the SK model in the presence of the 
above-mentioned random magnetic fields; a rich critical behavior is presented, and in par- 
ticular, one finds a critical frontier that may present one, or even two, tricritical points. 
The instabilities of the RS solution are also investigated, and AT lines presenting an in- 
fiection point, in concordance with those measured in some diluted antiferromagnets, are 
obtained. Finally, in section 4 we present our conclusions. 



II. THE ISING SPIN GLASS IN THE PRESENCE OF A RANDOM-FIELD 



The infinite-range-interaction Ising spin-glass model, in the presence of an external 
random magnetic field, may be defined in terms of the Hamiltonian 

H = —''^JijSiSj — '^^HiSi , (1) 

where the sum ^(jj) applies to all distinct pairs of spins Si = ±1 (z = 1,2, ...,N). The 
interactions {Jij} and the fields {Hi} follow independent probability distributions, 
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FIG. 1: The probability distribution of Eq. ([3]) (the random fields are scaled in units of a) for 
typical values of the ratio Ho/a: (a) (Ho/a) = 1/3,1,5/2; (b) (Ho/a) = 10. 
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The probability distribution for the fields {Hi} is a double Gaussian and depends on two 
parameters, i^o and a, modifying its form according to the ratio Hq/ct, as exhibited in 
Fig. 1. Such a distribution is double-peaked for (Ho/a) > 1, presents a single peak for 
{Hq/ct) < 1, changing its concavity at the origin when (Ho/a) = 1. Besides that, in the 
limit cr — > 0, one recovers a bimodal probability distribution. It is important to notice that 
its kurtosis, k =< Hf >h /[3(< H"^ >h)'^] [where < >h denotes an average over P(ifj)], 
varies from k = 1/3 (bimodal limit) to k = 1 (Gaussian limit), approaching unit only in 
the limit {Hq/cx) 0, in which case one gets a perfect Gaussian distribution. For finite 
values of H^/a one gets 1/3 < /t < 1, and in particular, for the cases exhibited in Fig. 1 
one has that k ^ 0.99 [(Ho/a) = 1/3], k ^ 0.83 [{Ho/ a) = 1], k ^ 0.50 [{Ho/a) = 5/2], 
and K ^ 0.35 [{Ho/a) = 10]. 

The free energy of a given disorder realization, F{{Jij, H^}), depends on two random 
variables, in such a way that the average over the disorder [ ]j^h niay be written in terms 
of independent integrals. 
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[F{{J,,,H,})]j^H = I \{[dJ,,P{J,,)]\{[dH,P{H,)]F{{.h,,H,}) . (4) 

Now, one can make use of the replica method [sl, I^, 0] in order to obtain the free energy 
per spin, 

-/3/= hm l[lnZ({J,„/J,})],,H= hm hm - 1) , (5) 

where represents the partition function of the rephcated system and f3 = l/{kT). 
Standart calculations lead to 

/3/ = - + lim imin g{m", , (6) 

4 2 n^o n 



where 



^7(m",g"^) = ^J2(^n' + - ^l^Tr.exp(7^+ ) - ^ In Tr„ exp(H;,) , 

(7) 



« (o/3) (a/3) « 



In the equations above, the index a {a = 1,2, ...,n) is a replica label, Tr^, represents a 
trace over the spin variables of each replica, and X](a/3) denote sums over distinct pairs of 
replicas. 

The extrema of the functional gim", give us the equilibrium equations 



= ^ < 5" >+ +^ < 5" >^ , (9) 
g"^ = i < >+ +1 < 5"^ >_ (ay^P) , (10) 
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where < >± indicate thermal averages with respect to the "effective Hamiltonians" 
in Eq. ^. 

Assuming the RS ansatz 

mm 

, i.e., = m (Va) and g"^ = q [V(a;/3)], Eqs. (E])- 

(fTOj) yield 
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where 



G 



q + 



1/2 



(14) 
(15) 



Although the spin-glass order parameter [Eq. ( |T3l) ] is always induced by the random 
field, it may still contribute to a nontrivial behavior. The RS solution is known to lead to 
an instability at low temperatures, usually associated to this parameter, occurring below 
the AT y line. 



(16) 



Let us now present the phase diagrams of this model. Since the random field induces 
the parameter q, there is no spontaneous spin-glass order, like the one found in the SK 
model. However, there is a phase transition related to the magnetization m, in such a way 
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that two phases are possible within the RS solution, namely, the Ferromagnetic (m 7^ 0, 
g 7^ 0) and the Independent (m = 0, g 7^ 0) ones. The critical frontier separating these 
two phases is obtained by solving the equihbrium conditions, Eqs. ( fT2l) and ( |T3i) . whereas 
in the case of first-order phase transitions, the free energy per spin, Eq. ffTTl) . will be 
analyzed. Expanding the magnetization [Eq. (fT2l) ] in power series, 

m = Ai{q)m + A3{q)m^ + A5{q)m^ + 0{m^) , (17) 



where 



A,{q) = /3Jo{l-pi(g)} , (18) 
AM = -^^{1 - 4pi(g) + 3p2{q)} , (19) 

A,iq) = i^{2 - 17pi(g) + 30p2(g) - 15p3(g)} , (20) 



and 



Pk{q) 



^ J^^ dze-'"/Hanh^' PJ !^Gz + . (21) 



The coefficients in Eqs. (fT8|) - (!20|) depend on q, which in its turn depends on m through 
Eq. f|T3|l . In order to eliminate this dependence, we expand Eq. f|T3|) in powers of m, 

q = qo + (f^'^o) \_^^j^,^ + O(m^) , (22) 



with 



r = l-4pi(go) + 3p2(go) , (23) 



where go corresponds to the solution of Eq. (|T3l) for m = 0. Substituting Eq. fl22l) in the 
expansion of Eq. (1171) . one obtains the m-independent coefficients in the power expansion 
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of the magnetization; in terms of the lowest-order coefficients, one gets, 



A[ = Ai{qo) , 



a: 



3 [ 1 - {(3jyT _ 



r . 



(24) 

(25) 
(26) 



The associated critical frontier is determined through the standard procedure, taking into 
account the spin-glass order parameter [Eq. ( |T3l) ]. as well. For continuous transitions, 
A\ = 1, with < 0, in such a way that one has to solve numerically the equation 
A'^ = 1, together with Eq. (fT3l) considering m = 0. If > 0, one may have first-order 
phase transitions, characterized by a discontinuity in the magnetization; in this case, the 
critical frontier is found through a Maxwell construction, i.e., by equating the free energies 
of the two phases, which should be solved numerically together with Eqs. ( fT2|) and ( |T3i) 
for each side of the critical line. When both types of phase transitions are present, the 
continuous and first-order critical frontiers meet at a tricritical point that defines the limit 
of validity of the series expansion. The location of such a point is determined by solving 
numerically equations A[ = 1, A'^ = 0, and Eq. (fT3ll with m = [provided that the 
coefficient of the next-order term in the expansion of Eq. ( 12^ is negative, i.e., < 0]. 

Considering the above-mentioned phases, the AT instability of Eq. (fT6l) splits each of 
them in two phases, in such a way that the phase diagram of this model may present four 



25 



phases, that are usually classified as 24 



I]: (i) Paramagnetic (P) (m = 0; stability 
of the RS solution); (ii) Spin-Glass (SG) (m = 0; instability of the RS solution); (iii) 
Ferromagnetic (F) (m ^ 0; stability of the RS solution); (iv) Mixed Ferromagnetic (F') 
(m 7^ 0; instability of the RS solution). 

Even though in most cases the AT line is computed numerically, for large values of 
Jo [i.e., Jq » J and Jq » Hq] and low temperatures, one gets the following analytic 
asymptotic behavior. 
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FIG. 2: Phase diagrams of the infinite-range-interaction ISG in the presence of a double-Gaussian 
random field; the phases are labelled according to the definitions in the text. ATI and AT2 
denote AT lines, and all variables are scaled in units of J. Two typical examples [(a) (c/ J) = 0.2; 
(b) (o"/J) = 0.6] are exhibited, for which there are single points (represented by black dots) 
characterized hy = 1 and ^3 = 0, defining the corresponding threshold values HQ^\a). 

For the particular case a = 0, i.e., the bimodal probability distribution for the fields 



25l | , it was verified that the phase diagrams of the model change qualitatively and quan- 
titatively for incresing values of Hq. Herein, we show that the phase diagrams of the 
present model change according to the parameters of the distribution of random fields 
[Eq. ([3])], which may modify drastically the critical line separating the regions with m = 
and m 7^ 0, defined by the coefficients in Eq. (^^. In particular, one finds numerically 
a threshold value, H^^\a), for which this line presents a single point characterized by 
A'i = 1 and A'-^ = 0; all other points of this line represent continuous phase transitions, 
characterized by A'^ = 1 and A'^ < 0. Typical examples of this case are exhibited in Fig. 2, 
for the dimensionless ratios (cr/J) = 0.2 and (cr/J) = 0.6. As will be seen in the next 
figures, for values of Hq/J slightly larger than H'^\a)/ J, this special point splits in two 
tricritical points, whereas for values of Hq/J smaller than H^\a)/ J, this critical frontier 
is completely continuous. Therefore, one may interpret the point for which Hq = ifg^'' (a) 
as a collapse of two tricritical points. Such an unusual critical point is a characteristic 
of some infinite-range-interaction spin-glasses in the presence of random magnetic fields 



25|, l26|], and to our knowledge, it has never been found in other magnetic models. /^From 
Fig. 2, one notices that the threshold value H'^\a)/J increases for increasing values of 
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(j/J, although the corresponding ratio HQ^\a)/(j decreases. Apart from that, this pecu- 
har critical point always occurs very close to the onset of RSB; indeed, for (cr/J) = 0.6, 
this point essentially coincides with the union of the two AT lines (ATI and AT2). At 
least for the range of ratios cr/J investigated, this point never appeared below the AT 
lines, i.e., in the region of RSB. Therefore, an analysis that takes into account RSB, will 
not modify the location of this point in these cases. 

In Fig. 3 we exhibit phase diagrams for a fixed value of o" (o" = 0.2 J), and increasing 
values of Hq. In Fig. 3(a) we show the case (Hq/J) = 0.5, where one sees a phase 
diagram that looks like, at least qualitatively, the one of the SK model; even though the 
random-field distribution [cf. Eq. ([3])] is double-peaked (notice that (Ho/a) = 2.5 in this 
case), the effects of such a field are not sufficient for a qualitative change in the phase 
diagram of the model. As we have shown above [see Fig. 2(a)], qualitative changes only 
occur in the corresponding phase diagram for a ratio {H^'^\a)/a) ^ 5, or higher. It is 
important to remark that a tricritical point occurs in the corresponding RFIM for any 



Ho/a) > 1 [30|, in agreement with former general analyses [31 



32 



If one associates 



the tricritical points that occur in the present model as reminiscents of the one in the 
RFIM, one notices that such effects appear attenuated in the present model due to the 



bond randomness, as predicted previously for short-range-interaction models [3J, |35|] . In 
Fig. 3(b) we present the phase diagram for {Hq/ J) = 0.993; in this case, one observes two 
finite-temperature tricritical points along the critical frontier that separates the regions 
with m = and m ^ 0. The higher-temperature point is located in the region where the 
RS approximation is stable, and so, it will not be affected by RSB effects; however, the 
lower-temperature tricritical point, found in the region of instability of the RS solution, 
may change under a RSB procedure. In Fig. 3(c) we exhibit another interesting situation 
of the phase diagram of this model, for which the lower-temperature tricritical point goes 
down to zero temperature, defining a second threshold value, Hq (cr). This threshold value 
was calculated analytically, through a zero-temperature approach that follows below, for 
arbitrary values of a/ J. Above such a threshold, only the higher-temperature tricritical 
point (located in the region of stability of the RS solution) exists; this is shown in Fig. 3(d), 
where one considers a typical situation with Hq > Hq (a). It is important to notice that 
in Fig. 3(d) the two AT lines clearly do not meet at the critical frontier that separates the 
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FIG. 3: Phase diagrams of the infinite-range-interaction ISG in the presence of a double-Gaussian 
random field with (cr/J) = 0.2 and typical values of Hq/J; the phases are labelled according 
to the definitions in the text. ATI and AT2 denote AT lines, and all variables are scaled in 
units of J. By increasing the value of Hq/J, the phase diagram changes both qualitatively 
and quantitatively and, particularly, the critical lines separating the regions with m = and 
m 7^ are modified; along these critical frontiers, the full (dotted) lines represent continuous 
(first-order) phase transitions and the black dots denote tricritical points; for the values of Hq/J 
chosen, one has: (a) continuous phase transitions; (b) two tricritical points at finite temperatures; 
(c) the lower tricritical point at zero temperature, defining the corresponding threshold value 

(2) 

Hq (a); (d) a single tricritical point at finite temperatures. 

regions with m = and m 7^ 0; such an effect is a consequence of the phase coexistence 
region, characteristic of first-order phase transitions, and has already been observed in 
the SK model with a bimodal random- field distribution The line ATI is valid up to 
the right end limit of the phase coexistence region, whereas AT2 remains valid up to the 
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FIG. 4: Phase diagrams of the infinite-range-interaction ISG in the presence of a double-Gaussian 
random field with (cr/J) = 0.6 and typical values of Hq/J; the phases are labelled according to 
the definitions in the text. ATI and AT2 denote AT lines, and all variables are scaled in units of 
J. Along the critical lines separating the regions with m = and m ^ 0, the full (dotted) lines 
represent continuous (first-order) phase transitions and the black dots denote tricritical points; 
for the values of Hq/J chosen, one has: (a) two tricritical points at finite temperatures; (b) the 
lower tricritical point at zero temperature, defining the corresponding threshold value Hq (<t). 

left end limit of such a region; as a consequence of this, the lines ATI and AT2 do not 
meet at the corresponding Independent-Ferromagnetic critical frontier. 

Additional phase diagrams are shown in Fig. 4, where we exhibit two typical cases for 
the random-field width (cr/J) = 0.6. In Fig. 4(a) we show the equivalent of Fig. 3(b), 
where two tricritical points appear at finite temperatures; now one gets qualitatively a 
similar effect, but with a random- field distribution characterized by a smaller ratio Ho/a. 
From the quantitative point of view, the following changes occur, in the critical frontier 
Independent-Ferromagnetic, due to an increase in a/ J: (i) such a critical frontier moves to 
higher values of Jq/ J, leading to an enlargement of the Independent phase [corresponding 
to the region occupied by the P and SG phases of Fig. 4(a)]; (ii) the two tricritical 
points are shifted to lower temperatures. In Fig. 4(b) we present the situation of a zero- 
temperature tricritical point, defining the corresponding threshold value Hq (ct); once 
again, one gets a physical situation similar to the one exhibited in Fig. 3(c), but with 
a much smaller ratio Hq/ct. Qualitatively similar effects may be also observed for other 
values of a, but with different threshold values, HQ^\a) and HQ^\a). We have noticed 
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FIG. 5: Evolution of the threshold values HQ^\a) (lower curve) and H^\a) (upper curve) with 
the width a (all variables are scaled in units of J). Three distinct regions (I, II, and III) are 
shown, concerning the existence of tricritical points and first-order phase transitions along the 
Independent-Ferromagnetic critical frontier. The dashed straight line corresponds to Hq = a, 



above which one has a tricritical point in the corresponding RFIM 



30t]. 



that such threshold values increase with cr/J, even though one requires less-pronounced 
double-peaked distributions [i.e., smaller values for the ratios Ho/a] in such a way to get 
significant changes in the standard SK model phase diagrams [as can be seen in Figs. 2, 
3(c), and 4(b)]. 

The evolution of the threshold values HQ^\a) and H^\(t) with the dimensionless 
width a/ J is exhibited in Fig. 5. One notices three distinct regions in what concerns 
the existence of tricritical points and first-order phase transitions along the Independent- 
Ferromagnetic critical frontier. Throughout region I [defined for Hq > if^^V)] a first- 
order phase transition occurs at finite temperatures and reaches the zero-temperature 
axis; a single tricritical point is found at finite temperatures [a typical example is shown 
in Fig. 3(d)]. In region II [defined for Hq^\o') < Hq < HQ^\a)] one finds two finite- 
temperature tricritical points, with a first-order line between them [typical examples are 
exhibited in Figs. 3(b) and 4(a)]. Along region III [Hq < HQ^\a)] one has a completely 
continuous Independent-Ferromagnetic critical frontier [like in Fig. 3(a)]. The dashed 
straight line corresponds to Hq = a, which represents the threshold for the existence of a 
tricritical point in the corresponding RFIM [30|]. Hence, if one associates the occurrence 
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FIG. 6: The zero-temperature phase diagram Hq versus Jq (in units of J) for two typical values of 
the dimensionless width o"/ J. The critical frontiers separating the phases SG and F' is continuous 
for small values of Hq/J (full lines) and become first-order for higher values of Hq/J (dotted 
lines); the black dots denote tricritical points. Although the two critical frontiers become very 
close near the tricritical points, they do not cross each other; the tricritical point located at a 
higher value of Jq/J corresponds to the higher dimensionless width a/ J. 

of tricritical points in the present model with those of the RFIM, one notices that such 
effects are attenuated due to the bond randomness, in agreement with Refs. 
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351; 



herein, the bond randomness introduces a spin-glass order parameter, in such a way that 
one needs stronger values of Hq/J for these tricritical points to occur. 

Let us now consider the phase diagram of the model at zero temperature; in this case, 
the spin-glass order parameter is trivial {q = 1), in such a way that the free energy and 
magnetization become. 
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where 



Gn 



1/2 



(30) 



Using a procedure similar to the one apphed for finite temperatures, one may expand 
Eq. (!29|) in powers of m, 



m = aim + asiri^ + a^mP + 0{m^) 



(31) 



where 




For [Hq/ {JGQy\^ < 1 [i.e., 03 < 0], we have a continuous critical frontier given by ai 



Jo /tT 



HI 



(35) 



This continuous critical frontier ends at a tricritical point (03 = 0), 



^f^l =1 



Ho 
J 



H, 



(2) 



J 



1/2 



(36) 



which may be substituted in Eq. ( l35l) to give 



Jo 
J 



lire 

Y 



1/2 



(37) 
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FIG. 7: Typical phase diagrams of tlie infinite-range-interaction ISG in tlic presence of a double- 
Gaussian random field with (cr/J) = 0.4 are compared with those already known for some 
particular cases. Comparisons of qualitatively similar phase diagrams are presented, essentially 
in what concerns the critical frontier that separates the regions with m = and m ^ 0. (a) 
Phase diagrams for the single Gaussian [{Hq/J) = 0.0 and (cr/J) = 0.4] and the double Gaussian 
[{Hq/J) = 0.8 and (cr/J) = 0.4] distributions for the random fields, (b) Phase diagrams for the 
bimodal [{Hq/J) = 0.9573] and the double Gaussian [{Hq/J) = 1.0447] distributions for the 
random fields, (c) Phase diagrams for the bimodal [{Hq/J) = 0.97] and the double Gaussian 
[{Hq/J) = 1.055] distributions for the random fields, (d) Phase diagrams for the bimodal 
[{Hq/J) = 1.0] and the double Gaussian [{Hq/J) = 1.077] distributions for the random fields. 
The phases are labelled according to the definitions in the text. ATI and AT2 denote AT lines, 
and all variables are scaled in units of J. 
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FIG. 8: Instabilities of the replica-symmetric solution of the infinite-range-interaction ISG (cases 
Jo = 0) in the presence of a double-Gaussian random field, for two typical values of distribution 
widths: (a) (cr/J) = 0.2; (b) {a/ J) = 0.6. In each case the AT line separates a region of RS 
from the one characterized by RSB (all variables are scaled in units of J). 

Hence, Eqs. ( 136!) and (137|) yield the coordinates of the tricritical point at zero temperature. 
In addition to that, the result of Eq. fl36|) corresponds to the exact threshold value H^^\a) 
(as exhibited in Fig. 5) . The above results are represented in the zero-temperature phase 
diagram shown in Fig. 6, where one finds a single critical frontier separating the phases 
SG and F'. 

In order to illustrate that the present model is capable of reproducing qualitatively the 
phase diagrams of previous works, namely, the Ising spin glass in the presence of random 



fields following either a Gaussian [2J|, or a bimodal [25|] probability distribution, in Fig. 7 
we compare typical results obtained for the Ising spin-glass model in the presence of a 
double Gaussian distribution characterized by (cr/J) = 0.4 with those already known for 
such particular cases. In these comparisons, we have chosen qualitatively similar phase 
diagrams, mainly taking into account the critical frontier that separates the regions with 
m = and m ^ 0. In Fig. 7(a) we exhibit the phase diagram of the present model 
[{Hq/ J) = 0.8] together with the one of an ISG in the presence of random fields described 
by a single Gaussian distribution; both phase diagrams are qualitatively similar to the one 
of the standard SK model. In Fig. 7(b) we present phase diagrams for the bimodal and 
double Gaussian distributions, at their corresponding threshold values, HQ^\a). Typical 
situations for the cases of the bimodal and double Gaussian distributions, where two 



tricritical points appear along the critical frontier that separates the regions with m = 
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and m 7^ 0, are shown in Fig. 7(c). Phase diagrams for the bimodal and double Gaussian 
distributions, at their corresponding threshold values, H'^\a)^ are presented in Fig. 7(d). 
Next, we analyze the AT instability for Jq = 0; in this case, Eq. (fT6l) may be written 

as 



(3i 



which corresponds to the same instability found in the case of a single-Gaussian random 
field I24I . In Fig. 8 we exhibit AT lines for two typical values of distribution widths; in each 
case the AT line separates a region of RS from the one characterized by RSB. One notices 
that the region associated with RSB gets reduced for increasing values of a; however, the 
most interesting aspect in these lines corresponds to an inflection point, which may be 
identified with the one that has been observed in the experimental equilibrium boundary 



of the compound FexZni_xF2 
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2J] . Up to now, this effect was believed to be explained 
only through the ISG in the presence of a single- Gaussian random field, for which the phase 
diagrams in the cases Jq > are much simpler, with all phase transitions being continuous, 
typically like those of the SK model. Herein, we have shown that an inflection point in 
the AT line may also occur in the present model, for which one has a wide variety of phase 
diagrams in the corresponding case Jo > 0, as exhibited above. Therefore, the present 
model would be appropriate for explaining a similar effect that may be also observed 
experimentally in diluted antiferromagnets characterized by first-order phase transitions, 
like FexMgi_xCl2. 



III. CONCLUSIONS 



We have studied an Ising spin-glass model, in the limit of infinite-range interactions 
and in the presence of random magnetic fields distributed according to a double-Gaussian 
probability distribution. Such a distribution contains, as particular limits, both the single- 
Gaussian and bimodal probability distributions. By varying the parameters of this distri- 
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bution, a rich variety of phase diagrams is obtained, with continuous and first-order phase 
transitions, as well as tricritical points. The condition for the existence of a single finite- 
temperature tricritical point at the critical frontier Paramagnetic- Ferromagnetic (i.e., in 
the region of stability of the replica-symmetric solution), characterized by a first-order line 
at low temperatures, is derived analytically. Besides that, we found an infiection point in 
the AT line (in the plane magnetic field versus temperature), which may correspond to the 



one observed in the compound FexZni_xF2 
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24| . This effect, which has a. 



ready shown 



up in the Ising spin-glass model in the presence of a Gaussian random field [2J], is herein 
obtained for a more general probability distribution for the magnetic fields. Hence, the 
present model is appropriate for explaining a similar effect that could be observed also in 
diluted antiferromagnets characterized by first-order phase transitions, like FexMgi_xCl2. 
Therefore, with this random-field distribution, one may adjust the model to given physical 
situations, in order to reproduce a wide diversity of effects that occur in real systems. 

The double-Gaussian probability distribution, defined above, is suitable for a theo- 
retical description of random-field systems, being a better candidate for such a purpose 
than the two most commonly used distributions in the literature, namely, the bimodal and 
single-Gaussian distributions, due to the following reasons: (i) In the identifications of the 
RFIM with diluted antiferromagnets in the presence of a uniform magnetic field, the local 
random fields are expressed in terms of quantities that vary in both sign and magnitude 
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291]. This characteristic rules out the bimodal probability distribution from such a 



class of physical systems. The double-Gaussian probability distribution is appropriate for 
a description of diluted antiferromagnets for a large range of magnetic concentrations, like 
in the RFIM, as well as in the ISG regimes, (ii) Although the RFIM defined in terms of a 
single-Gaussian probability distribution for the fields is physically acceptable, it usually 
eads to a continuous phase transition at finite temperatures, either within mean-field 
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331], or standard short-range-interaction approaches (36|, l37|. Such a system is not 



able to exhibit first-order phase transitions and tricritical points, that may occur in some 



diluted antiferromagnets 



141 ]. A similar behavior was obtained for an ISG in the presence 



of random magnetic fields following a single- Gaussian probability distribution, where all 
phase transitions were found to be continuous 2J]. The present model, defined in terms 
of a double- Gaussian probability distribution, is expected to be relevant for FexMgi_xCl2 
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(which is known to exhibit a first-order phase transition in the RFIM regime 13]) with 
concentrations x < 0.55, where the ISG behavior shows up. 
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